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ABSTRACT 

The effect of a dynamo-generated mean magnetic field of Beltrami type on the mean electromotive force is 
studied. In the absence of the mean magnetic field the turbulence is assumed to be homogeneous and isotropic, 
but it becomes inhomogeneous and anisotropic with this field. Using the testfield method the dependence of the 
alpha and turbulent diffusivity tensors on the magnetic Reynolds number Re^ is determined for magnetic fields 
that have reached approximate equipartition with the velocity field. The tensor components are characterized 
by a pseudoscalar a and a scalar turbulent magnetic diffusivity rj,. Increasing Re^ from 2 to 600 reduces rj, by a 
factor X! 5, suggesting that the quenching of 77, is, in contrast to the 2-dimensional case, only weakly dependent 
on Rcjvf. Over the same range of Re^, however, a is reduced by a factor « 14, which can qualitatively be 
explained by a corresponding increase of a magnetic contribution to the a effect with opposite sign. The level 
of fluctuations of a and 77? is only 10% and 20% of the respective kinematic reference values. 
Subject headings: MHD - turbulence 
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1. INTRODUCTION 

Magnetic fields in stars and galaxies tend to display large 
scale spatial order, and in the case of the Sun also long term 
temporal order (the 22 year cycle). The underlying process 
is generally believed to be a turbulent large-scale or mean- 
field dynamo - the simplest of which is an dynamo, which 
works with helical turbulence and no mean flows. This can 
be modeled by direct numerical simulations in a periodic box 
where the flow is driven by helical isotropic forcing. Corre- 
sponding simulations by Brandenburg (2001) show that in the 
nonlinear regime there is a resistively slow saturation phase 
associated with nearly perfect conservation of magnetic he- 
licity. This slow saturation imposes tight constraints on the 
quenching of the electromotive force. By comparing with 
suitable mean field models one can only constrain the quench- 
ing of the full electromotive force, but not the individual 
quenchings of a and 77,, because the saturated mean mag- 
netic field of an a~ dynamo tends to become force-free, so the 
mean magnetic field and the mean current density are aligned 
(Blackman & Brandenburg 2002; hereafter BB02). As a con- 
sequence an infinitude of combinations of quenching expres- 
sions for a and rj, describe the same saturation behavior. 

The saturation of the mean magnetic field is well described 
by a mutual cancellation of kinetic and magnetic alpha ef- 
fects, where the latter depends on the production rate of mean 
magnetic helicity. To reproduce the resistively slow satura- 
tion, both kinetic alpha effect, ok, and turbulent magnetic 
diffusivity, 77,, could be assumed completely unquenched. 
This is however an unrealistic simplification (Kleeorin & Ro- 
gachevskii 1999). Some level of quenching of 77, was found 
to be necessary to reproduce the simulations (BB02). 

Since the early work of Vainshtein & Cattaneo (1992), a 
lot of effort has gone into determining the quenching of a. 
It is now clear that for mean fields defined as volume aver- 
ages over a periodic box a is "catastrophically" quenched like 
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Re^ with mean fields of equipartition strength (Cattaneo & 
Hughes 1996). However, subsequent work showed that this is 
a particular consequence of the use of full volume averages, 
in which case the mean current density is zero (BB02). 

The quenching of 77, is much less understood. While 
in the two-dimensional case, 77, is indeed catastrophically 
quenched (Cattaneo & Vainshtein 1991), in three dimensions 
the quenching may depend just on B^, but not on Re^. This 
has already been found from the decay rate of a nonheli- 
cal large-scale magnetic field in driven non-helical turbulence 
(Yousef et al. 2003). Similar indications come also from fit- 
ting mean field models to corresponding simulations (BB02). 

Quantifying more precisely the simultaneous quenching of 
a and 77, is the goal of the present paper We admit both a and 
77t to be tensors, denoted by ajj and 77,^, respectively, and we 
calculate them using the testfield method (e.g., Brandenburg 
et al. 2008, Suret al. 2008). However, unlike earlier kinematic 
work, we now allow the velocity to be the result of the fully 
nonlinear hydromagnetic equations, i.e. to be influenced by 
the resulting mean magnetic field. 

2. THE METHOD 

Following earlier work by Brandenburg (2001), we con- 
sider a compressible isothermal gas with sound speed c,, but 
in addition we also solve a set of testfield equations, as was 
done in Brandenburg et al. (2008) for the kinematic case. The 
full set of governing equations is then 



^ = -U-VU-c^-^Xnp + f+p-^i^JxB + V-lpvS^, 
d\np 



dt 



= -{7-Vlnp-V-{7, 



dA_ 
~dt 



UxB - jU()T]J, 



(1) 



(2) 



(3) 



daPi - 



dt 



^UxbP'' + uxBi"^ + uxb'"' -ux bi'i - i^orif, (4) 



where mean fields are defined as horizontal (xy) averages, 
thus being functions of z and t only, and indicated by over- 
bars whereas lower case vectors denote deviations from the 
averages ("fluctuations"). The superscripts pq refer to four 
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separate equations that are characterized by four different test- 
fields BP'' having a cos kz or sin kz dependence {q - c, s) in 
the jc or y component {p - 1,2). We employ a magnetic vec- 
tor potential both for the magnetic field B - V x A and for 
the responses to the testfields, = V x a'"'. We reinitial- 
ize a''^ to zero every 30-60 turnover times to suppress small- 
scale dynamo action (cf. Sur et al. 2008). Of course, the ve- 
locity U is now affected by the magnetic field B through the 
Lorentz force. The current density is / = V x Bj^Q, where yUo 
is the magnetic permeability. The flow is driven by random 
forcing described by a forcing function / consisting of circu- 
larly polarized plane waves with positive helicity and random 
direction (giving rise to a flow with maximal helicity), and 
Sij - jiUij + Ujj) - l^/jV ■ U is the traceless rate of strain 
tensor The forcing function is chosen such that the moduli 
of the wavevectors, \kf\, are in a narrow interval around an 
average value, which is denoted simply by kf. 

Owing to our definition of averages, B is independent of x 
and y and all its first-order spatial derivatives can be expressed 
by the components of 7. If we ignore higher-order derivatives 
of B the mean electromotive force £i - u xb has the form 

6i^aijBj-iUQT]ijlj (5) 

with two tensors aij and 77,^, and we restrict our attention 
to 1 < /, j < 2. For details see Brandenburg et al. (2008). 
Solving the test field equations allows us to calculate - 
u X bPt and, via Eq. (|5]), all 4+4 components of aij and rfij. 

Important control parameters are the magnetic Reynolds 
and Prandtl numbers, Rcm - u^m^lirjkf) andPr^ - v/?;, where 
Mims = {u^y^^ is the actual (magnetically affected) rms veloc- 
ity and angular brackets denote volume averages. The small- 
est possible wavenumber in a triply-periodic domain of size 
L X L X L is ki = 2n/L. In order to achieve large values of 
Re^f, the value of kf/ki should be small, but still large enough 
to allow for a clear separation of scales between the domain 
scale and the energy-carrying scale. We use kf/ki = 3 as a 
compromise. 

The structure of the turbulence is determined by the vectors 
B and /, but for a Beltrami field they are aligned, so we have 

aij(B) = ai (B)6ij + a2(B)BiBj, (6) 

7],j(B) = m(B)6ij + m(B)BiBj, (7) 

where B means the unit vector in the direction of B. When 
inserting this into the general expression for the electromo- 
tive force given above this reduces to £ = aB - ^oT],] , with 
coefficients 

a - a\ + a2 - ri2kn and rj, - rju (8) 

where k^ - k„{z, t) = hqJ ■ B/B^ isa pseudoscalar that quan- 
tifies the helicity of the large-scale field. (Here ki„/ki x -I.) 
We emphasize that for Beltrami fields the assignment of ai , 
a2, r]\ and 772 to a and 771 is not unique. In the general situa- 
tion, when the mean field is not of Beltrami type, instead of a2 
and 772 eight new coefficients emerge which contribute in an 
unambiguous way to field generation and dissipation. Future 
work must show whether our a\ and 77, are then still dominant. 

3. RESULTS 

Throughout this paper we fix Pr^ = 1 and vary Re^ be- 
tween 2 and 600. For large values of Re^ a broader range 
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Fig. 1 . — Compensated time-averaged spectra of kinetic and magnetic 
energy, as well as of kinetic and magnetic helicity. for a run with Re^ = 600. 
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Fig. 2. — Visualization of B, on the periphery of the compu- 
tational domain for a run with Re^ = 600 and a resolution of 
512' mesh points. Note that on average the field is compatible 
with that in equation 0. Note also the clear anisotropy with struc- 
tures elongated in the direction of the field. For an animation see 
http: //www.nordita.org/so£tware/pencil-code/movies/icascade/| 

of scales is excited, as can be seen in spectra of kinetic and 
magnetic energy, Ek and Em, shown in Fig. [1] In the range 
4 < k/ki < 30 both spectra are comparable to a k^^^^ spec- 
trum. For comparison, spectra of kinetic and magnetic helic- 
ity, Hk and Hm, are also shown. 

For Rcm < 2 there is no dynamo action, but in all other 
cases a large scale magnetic field is maintained (Fig. |2]i, just 
as in Brandenburg (2001), except that here kf/ki - 3 instead 
of 5 or larger The dynamo is of type and hence the mean 
field a Beltrami field, 

B(z,t) ^B(t){cos0,sm0,O), d^kiz + (p, (9) 

with phase <p. To shorten the transient phase we use this field 
also as initial condition. 

Inserting (|9]l into (|6]l and (|7]i and calculating suitable aver- 
ages over z (or volume) we get 

a2 = 8 {ai2 cos 0sin0) = 8 (q'21 cos ^?sin0) , (10) 
ai + a2/2 = (cii) = <a22) , (H) 
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and analogous for rji and 772. Obviously, the determination of 
ai, ai, rii, and 7/2 requires knowledge of the Beltrami phase (p, 
which often drifts away from its initial value during the course 
of the run. We determined therefore the actual phase 0(f) by 
applying a suitable Fourier analysis to B. 

In general, a quenching can involve time derivatives (e.g., 
Kleeorin & Ruzmaikin 1982, BB02). In order to avoid such 
complications we focus on statistically steady (dynamo) solu- 
tions, that is, on the saturated dynamo fields. For given values 
of the parameters of the system ([T])-(l3|, the saturation strength 
of B is uniquely determined. Hence, by changing the forcing 
strength or 77 we are only able to follow a specific path in the 
B - Rcm plane, but not to scan it in a 2D fashion. 

In Table [1] we represent the results in nondimensional form 
with normalized quantities indicated by a tilde. We normalize 
the rms values of the mean field and the fluctuations with the 
equipartition field strength Z?eq - (/^o(pM^(^)))'^^ and intro- 
duce 

m^mhto^ m^m/^to, fj^ri/m' (12) 

&i - ai/ao, 02 - ailoQ, (13) 



1.00 



where 77,0 = \u,^J,B)k^.^ , 



3^rm,sv"/-/ , Co = -iMrms(B), and u,^JB) is 
the rms velocity of the saturated state, so the reference values 
are already magnetically affected. This normalization implies 
that in the kinematic case ai = ^1 = 1 (Sur et al. 2008), 
while Q'2 - fji - 0. Error bars are calculated based on the 
maximum departure obtained from the three time series, each 
taken over one third of the full sequence. 

The consistency of the results for a and with the presence 
of a steady state can be assessed by calculating the growth 
rate, A, of the associated kinematic mean field dynamo for a 
Beltrami field with k,,, - -k\, i.e. A - -ak\ - (rjt + T])ky 
In the saturated state A should vanish. Again, we present A 
in nondimensional form, here in terms of the turbulent decay 
rate, 



A = Aliriiok^) = akf - {i), + fj). 



(14) 



where a - a/ao and kf = kf/ki. Within error bars, the value 
of A is consistent with zero, thus supporting the consistency 
of a and rj, with the established steady state; see Table[T] (An 
exception is Run A, because it is subcritical and so /I < 0.) 
This in turn supports the applicability of the testfield method 
to the nonlinear case. However, as in almost all supercritical 
runs a small-scale dynamo is operative, our results which are 
derived under the assumption of its influence being negligible 
may contain a systematic error If present, it should be small 
though, given the good precision of the results for A. A more 
thorough study of the role of the small-scale dynamo will be 
the subject of future work. 

A measure of the reliability of the averages is the length of 
the time series in "turnover" times. At = Mi-ms^f(fmax - fmin)- 
Our results presented in Table [T] show a decline of a by a fac- 
tor » 15 and a decline of 77, by a factor w 5 as Rcm increases 

by a factor 300 while B = 0.95...1.4Beq. 

As expected, there are random fluctuations of a and rj,, rep- 
resented here by their non-dimensional rms values, o-rms - 
Q-rms/ao and 7],^^ = VimJlto- Even for large Re^ the fluctu- 
ations remain around 0.1 and 0.2, respectively. This is less 
than in the kinematic case (Brandenburg et al. 2008), but still 
comparable to the mean values of a and 77/, respectively. 

4. DISCUSSION 
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Fig. 3. — ReM-dependence of a and fji/kf together with and -S*/. 

Let us now put our results in relation to earlier work, which 
mostly used mean fields defined as full volume averages, 
hence being uniform. In that case a was quenched all the 
way to zero like ReJi^' . This result can be understood in terms 
of a mutual cancelation of kinetic and magnetic contributions 
to the a effect (Pouquet et al. 1976), 

a ^aK + om, Ok = ■ «, qm = 5TJ • fe/p, (15) 

where (o - V x «. Assuming Tu^^skf ~ 1 (Brandenburg & 
Subramanian 2007), we estimate r and hence, by measuring 
(co-u) and (J-b), we determine a/c = ag/a^) and um = aM/o;^', 
see Table[T]and Fig. [3] It turns out that uk is essentially inde- 
pendent of Rcm [but of course dependent on B; see Table 1 of 
Brandenburg (2001)] and om approaches a certain fraction of 
aiK, reducing the residual a in equation ( fTSb as Re^ increases. 
This agrees only qualitatively with the measured decline of a, 
because the residual a is sill too big. However, Eq. (fTSb as- 
sumes isotropy and that the values of t are the same for 
and om, which is not borne out by simulations (Brandenburg 
& Subramanian 2007). By contrast, our direct calculations 
show that a is quenched to values of order fji/kj , as is neces- 
sary for a steady state; see Eq. ( fT4b . Note that the decline of 
fit is much weaker than in the two-dimensional case where 77, 
decreases like Re^' (Cattaneo & Vainshtein 1991). 

5. CONCLUSIONS 

For the first time it has been possible to determine both aij 
and rjij in the magnetically quenched case. These tensors are 
here characterized by the non-tensorial quantities a and rji. 
The consistency of the results of the testfield method suggests 
that the nonlinear a can be determined by the knowledge of 
just u{x, t) over the past several correlation times-even if it is 
already influenced by B. Qualitatively, the quenching of a can 
be explained by -um approaching ok for finite field strengths 
and large Re^. Generally, a will be quenched to whatever is 
the value of (77, + Tf)k\ (BB02). However, until now we had 
no idea how big the quenched value of 77, is. There was the 
possibility that 77, was quenched to very small values, just like 
in the two-dimensional case (Cattaneo & Vainshtein 1991). If 
that were true, a would also be very small. We can now say 
that this is not the case, because 77, is only reduced to about 
20% of the kinematic value, while the normalized value a is 
quenched to ^ fi,/kf » 7% of its kinematic value, as is seen in 
Fig.S 
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TABLE 1 

Transport coEFPicmNTS for runs in the ranoe 2 < Rbm ^ 600 at saturation field strengths. 



Run 


ReAf 






a 


7< 




i 


-012 


-m 


Q'rms 


'^rms 




-Sim 


A? 


A 


2 


0.0 


0.0 


0.70 ± 0.03 


0.67 ± 0.07 


1.57 


-0.14 ±0.01 


0.04 ± 0.05 


-0.02 ± 0.06 


0.09 


0.12 


1.03 


0.01 


150 


B 


4 


0.9 


0.4 


0.44 ± 0.01 


0.58 ± 0.04 


0.73 


0.00 ± 0.00 


0.33 ± 0.02 


-0.11 ±0.03 


0.10 


0.21 


1.02 


0.31 


422 


C 


12 


1.7 


0.7 


0.24 ± 0.01 


0.46 ± 0.02 


0.25 


0.00 ± 0.00 


0.37 ± 0.02 


-0.04 ±0.01 


0.09 


0.16 


1.00 


0.55 


601 


D 


30 


1.9 


0.8 


0.16 ±0.01 


0.36 ± 0.02 


0.11 


-0.00 ± 0.01 


0.37 ± 0.02 


0.03 ± 0.03 


0.07 


0.14 


1.02 


0.62 


350 


E 


60 


2.0 


0.8 


0.09 ± 0.01 


0.22 ± 0.02 


0.05 


0.00 ± 0.01 


0.33 ±0.01 


0.05 ± 0.01 


0.09 


0.22 


1.00 


0.66 


711 


F 


150 


2.0 


0.9 


0.07 ± 0.00 


0.19 ±0.01 


0.02 


0.01 ± 0.01 


0.24 ± 0.05 


0.08 ± 0.01 


0.07 


0.16 


1.01 


0.69 


225 


G 


300 


1.8 


0.9 


0.06 ± 0.00 


0.15 ±0.00 


0.01 


0.01 ± 0.01 


0.21 ± 0.02 


0.05 ± 0.02 


0.06 


0.16 


1.01 


0.66 


177 


H 


600 


1.8 


0.9 


0.05 ± 0.01 


0.13 ±0.01 


0.005 


0.01 ± 0.04 


0.14 ±0.05 


0.04 ± 0.01 


0.05 


0.10 


1.03 


0.64 


44 



Obvious extensions of this work include the application to 
non-Beltrami fields and to domains with boundaries and/or 
shear. In the latter case there exists a great deal of earlier work 
with relevant simulation data supporting the idea of an a effect 
that is strongly controlled by magnetic heUcity evolution, and 
that catastrophic quenching can be decisively alleviated in the 
presence of shear-driven magnetic heUcity fluxes. 
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